Polarization dependence and symmetry analysis in indirect ^-edge RIXS 
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We present a study of the charge-transfer excitations in undoped Nd2Cu04 using resonant inelastic X-ray 
scattering (RIXS) at the Cu A^-edge. At the Brillouin zone center, azimuthal scans that rotate the incident- 
photon polarization within the CUO2 planes reveal weak fourfold oscillations. A comparison of spectra taken in 
different Brillouin zones reveals a spectral weight decrease at high energy loss from forward- to back-scattering. 
We show that these are scattered-photon polarization effects related to the properties of the observed electronic 
excitations. Each of the two effects constitutes about 10% of the inelastic signal while the '4p-as-spectator' 
approximation describes the remaining 80%. Raman selection rules can accurately model our data, and we 
conclude that the observed polarization-dependent RIXS features correspond to Eg and Big charge-transfer 
excitations to non-bonding oxygen 2p bands, above 2.5 eV energy-loss, and to an £g rf — > rf excitation at 
1.65 eV. 
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Raman scattering and optical spectroscopy have enabled 
tremendous contributions to the study of condensed matter 
systems. Both probes use ~ 1 eV Ught and are Hmited to es- 
sentially zero momentum transfer. In order to investigate the 
charge response of a material throughout the Brillouin zone, 
photons in the X-ray regime must be used. X-ray Raman scat- 
tering, more commonly referred to as resonant inelastic X-ray 
scattering (RIXS), allows the measurement of the momentum 
dependence of charge excitations. Even though it has suc- 
cessfully been used to study the physics of a wide array of 
systems,'"^ this resonant technique is relatively new and there 
is still much to learn about the details of its cross section. 

One of the strengths of conventional Raman scattering de- 
rives from the fact that the technique allows the determination 
of the symmetry of excitations by selecting the polarization 
of the incident and scattered photons. Only polarization de- 
pendent soft (or direct'') X-ray RIXS studies, at the O K-edge 
or the Cu L- and M-edges of the cuprates,— — have shown 
that photon polarization can be used to select different elec- 
tronic excitations. In these cases, the photon polarization ef- 
fects are understood to come from the direct interaction of the 
photoelectron, excited into the valence system, with the va- 
lence electrons. This is why this type of RIXS is commonly 
refeiTed to as direct.*" On the other hand, tuning the incident 
energy to the Cu K-edge excites the photoelectron into the Cu 
4p band approximately 10-20 eV above the 3d valence lev- 
els. The relatively large separation in energy from the valence 
levels as well as the large spatial extent of the 4p orbital are 
believed to prevent the photoelectron from interacting with 
the valence system, and this type of RIXS is in turn referred 
to as indirect.* Accordingly, theoretical models of this type 
of RIXS accordingly do not include possible interactions of 
the 4p photoelectron with the valence system, '^"'^ which is 
commonly referred to as the '4/7-as-spectator' approximation. 

RIXS investigations of the photon polarization effects have 
so far concluded that the incident polarization does not af- 
fect the probed valence band excitations, but only determines 
their specific resonance energies based on the crystal field lev- 
els of the 4p photoelectronr^ii^iii in accordance with the 4p- 



as-spectator approximation. Nonetheless, there still exists no 
quantitative experimental evidence for an incident-photon po- 
larization independence. Furthermore, the scattered-photon 
polarization effects at the transition metal K-edge have re- 
ceived very little attention, and it remains unclear whether 
they are present and can be used to learn about the symme- 
try of charge-transfer excitations in transition metal oxides. 

This paper is separated in three parts. In Section U we 
discuss the photon-polarization dependence expected within 
the 4/?-as-spectator approximation. In Sec. [Ill we study the 
photon-polarization and scattering-geometry dependence of 
RIXS at the Cu K-edge of Nd2Cu04 and observe scattered- 
photon polarization effects beyond the 4/:i-as-spectator ap- 
proximation. In Sec. Ill Al we present data obtained upon rotat- 
ing the incident-photon polarization within the CUO2 planes 
while in Sec. Ill B I we present a comparison of zone-center 
spectra taken in different Brillouin zones. The results are 
discussed in Sec. |lll] The normalization procedure used to 
correct for sample self-absorption and compare RIXS sig- 
nal across different scattering geometries is described in Ap- 
pendixlAl 



I. POLARIZATION DEPENDENCE OF THE RIXS CROSS 
SECTION 

Before presenting our experimental observations of photon 
polarization effects in RIXS, it is instructive to discuss what 
is expected according the current understanding of the indi- 
rect RIXS cross section. At the Cu A'-edge, the RIXS pro- 
cess starts with the resonant absorption of an X-ray photon 
which creates a I5 core hole and a 4p photoelectron on the 
Cu site. During the lifetime of the core hole, before it recom- 
bines with the photoelectron and an X-ray photon is emitted, 
the core-hole interacts strongly with the valence system to cre- 
ate electronic excitations. The 4p photoelectron, on the other 
hand, is believed to be only a spectator and to evolve without 
interacting with the valence system. 

In order to determine the cross-sectional implications of 
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this 4/?-as-spectator approximation, we consider a simple 
model where the photoelectron's energy eigenstates are dic- 
tated by the crystal-field: £"J!^ = A and eV' = 0. While this 
approach neglects Ap band eff'ects, in practice, the RIXS in- 
tensity is often confined within narrow and well-separated in- 
tervals in incident energy which should be approximated well 
by this two-level model. Since the initial and final states of 
the Ap photoelectron are determined by the initial (e,) and 
scattered {ef) photon polarizations, the transition amplitude 
(/4p) of the photoelectron is determined by its evolution within 
the crystal-field during the lifetime of the core hole (/i/F) and 
contains polarization dependent factors that modulate the in- 
tensity of the excitations created by the core-hole: 
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In the case where A » F, the cross section simplifies to a 
polarizer: only the polarization components along the res- 
onating 4/7 crystal-field level contribute to the inelastic signal: 



for an in-plane resonance (£, = A) and le^ef i 



for an out-of-plane resonance (£, = 0). 

In their study of CuO, Doring et al.^ pointed out that, 
within the 4/9-as-spectator approximation, the electric dipole 
absorption-emission matrix element is equivalent to the reso- 
nant elastic X-ray scattering cross section described in detail 
by Hannon et al. Furthermore, they successfully apply this 
approximation to describe the scattering angle dependence of 
the intensity of the 5.4 eV local charge-transfer excitation in 
horizontal scattering geometry. Equation[T]is a simplification 
of the formula presented in Ref. but it succinctly captures the 
important effect of scattering geometry on the RIXS signal. 
Including diff'erent valence states, like the well- and poorly- 
screened intermediate states, would improve it. 

In general, Eq. [1] shows that the inelastic intensity is max- 
imized when both the incident- and scattered-photon polar- 
izations are parallel to each other and point along a crystal- 
field eigenstate. These conditions can naturally be fulfilled in 
vertical-scattering geometry where both polarizations are per- 
pendicular to the scattering plane (cr polarized),^-'' For hori- 
zontal scattering geometry, with both polarizations within the 
scattering plane (n polarized), these conditions can only be 
approached for forward and backward scattering. 

On the other hand, in horizontal-scattering geometry the 
elastic intensity can be minimized independently of the in- 
elastic intensity. This can be an advantage since it allows 
the suppression of the elastic 'tail' due to the non-zero en- 
ergy resolution, leading to a higher signal-to-background ra- 
tio. At a scattering angle (26) of 90°, the non-resonant elastic 
contribution is zero because the polarization factor of Thomp- 
son scattering (ey ■ e,) is zero. To have non-zero inelastic 
intensity, the non-degenerate 4p eigenstates can be used as 
cross-polarizers. By polarizing the incident photon between 
4p crystal-field eigenstates of different energies, e.g., e, || jch-z, 
the inelastic signal will be dominated by the excitations cre- 



ated by the incident photon polarization component along the 
resonating 4p crystal-field eigenstate. Since this resonating 
component is not perpendicular to the scattered-photon polar- 
ization, the resonant inelastic signal will be detectable while 
the non-resonant elastic signal will be zero. Note that this ef- 
fectively rotates the photon polarization 90°. This should not 
be mistaken for a scattered-photon polarization effect where 
the 4p photoelectron scatters during the RIXS process. 

This type of scattering geometry is used extensively.i2Si21 
Based solely on the 4p polarization eff'ects, the maximum 
RIXS intensity in this horizontal scattering geometry is ob- 
served when both incident and scattered photon polarization 
are at 45° from the resonantly excited 4p eigenstate. In 
the tetragonal crystal-field symmetry of cuprates for exam- 
ple, 211^ ™d 26 - 90° allows maximum RIXS intensity for 
both in-plane and out-of-plane resonances while minimizing 
the non-resonant contribution to the elastic line. As is shown 
in Appendix lAl self-absorption effects will move this maxi- 
mum of RIXS intensity towards a more grazing incidence an- 
gle (keeping 26 - 90°) when the scattering surface is perpen- 
dicular to the [0 L] direction. 



II. PHOTON POLARIZATION EFFECTS IN NCO 

For our study of the scattering-geometry and photon- 
polarization dependence of RIXS, we chose Nd2Cu04 
(NCO), the tetragonal Mott-insulating parent compound 
of the electron-doped high-temperature superconductor 
Nd2-.vCevCu04. A single crystal was prepared as described 
previously^'* and studied in its as-grown state. A larger piece 
from the same growth was measured with neutron scattering 
and found to have a Neel temperature Tff ~ 270 

The RIXS data were collected in vertical scattering geom- 
etry with the X-ray spectrometer on beamline 9-ID-B, at the 
Advanced Photon Source, and using the 2 m arm configura- 
tion. The energy resolution (FWHM) was 0.32 eV (Sec. Ill Al l 
and 0.25 eV (Sec. HTBT i. 

The 'elastic tails' due to the elastic peak and non-zero en- 
ergy resolution are subtracted from the inelastic spectra above 
1.5 eV energy loss. This is accomplished by (i) using the elas- 
tic peak to establish zero energy transfer, (ii) fitting energy- 
gain data (typically up to 3 eV) to the heuristic modified 
Lorentzian form 1/(1 + Iwf) where w is energy-loss (with 
typical values of a in the range 1-2), (iii) and then subtracting 
the result of the fit from the energy-loss part of the spectrum. 



A. Azimuthal (i//-) scans 

In order to study the effects of photon polarization inde- 
pendently of the momentum transfer, we perform azimuthal 
scans, which rotate the incident-photon polarization within 
the Cu02 planes, while keeping the energy- and momentum- 
transfer constant (the rotation axis is parallel to Q.^^) This 
scattering geometry is illustrated in Fig. [T^, with at - 0, and 
in Fig. [TJ). The azimuthal angle i/r is the scanned variable. 
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FIG. 1. (a) Side and (b) top view of the scattering geometry used 
thiroughout tliis paper. The polarization of the incident beam (cr) and 
the two possible polarization conditions for the scattered beam (cr 
and n) are emphasized in both panels. Note that (/'-scans are collected 
at ai = only (b) inset: two-dimensional (2D) Brillouin zone: the red 
line shows the full region of integration arising from the momentum 
resolution whereas the blue line shows its half width, (c) Energy 
scans taken at the 2D zone center with incident energy E,=&997 eV 
(circles) and 9001 eV (squares). The superposed colored circles and 
square show where the tf/scms of corresponding color and energy 
are measured in (d) are measured. The charge-transfer gap mea- 
sured with optical conductivity is indicated for comparison.-* (d) 
Azimuthal scans taken in a full circle around the 2D zone center 
at 1.625, 2, 3.5, 5 eV energy-loss (for E, = 8997 eV) and at 6 eV 
energy-loss for (£i=9001 eV) for Q = (0 9.1) and symmetrized 
with respect to the 90° rotations and mirror planes of the underlying 
tetragonal structure, (e) Energy scans with elastic tail subtracted (see 
text) taken at g = (0 0.5 11.1) for different values of the azimuthal 
angle (tf/). The scans indicated by circles were taken at Ei = 8997 eV 
whereas those indicated by squares were taken at E^ = 9001 eV. En- 
ergy scans with tf/ = 45° are corrected for changes in self-absorption 
(+5.2%) and incident polarization (+0.8%) compared to = 10.7°. 
(f) Fourier components of the i^-scan at 3.5 eV energy-loss with am- 
plitude given in percent of 0-fold (DC) component (not shown). The 
4-fold component is well above the statistical noise level (gray re- 
gion). There is also a 2-fold component present in the raw data. 



Figure[T}; presents characteristic RIXS spectra taken at two 
different incident energies. These line scans show electronic 
excitations in the 1-10 eV energy-loss range associated with 
the electronic structure of the strongly-correlated Cu02 plane, 
and they indicate the energies where i/'-scans are made. At the 
incident energy Ei = 8997 eV, the chosen values span from 



1 .625 eV, just above the optical charge-transfer gap,^^ to the 
maximum of the inelastic intensity at 5 eV. At Ei = 9001 eV, 
we study the molecular orbital excitation^l by measuring at 
6 eV. The raw azimuthal scans are symmetrized to be con- 
sistent with the underlying tetragonal symmetry of the CuOa 
planes. This symmetrization consists of folding back the data 
within the 0-45° arc and then averaging it, the results are 
shown in Fig. [TJl. This procedure averages out all non-4-fold 
components exactly which, in practical terms, filters out ex- 
perimental noise and leaves only the intrinsic electronic prop- 
erties. For Ei - 8997 eV, the resulting amplitudes of the 4- 
fold oscillations ai-e -0.8 + 3.5% (at 1.625 eV), -0.2 + 0.9% 
(2 eV), 3.2 + 0.9% (3.5 eV), and 3.2 + 0.9% (5 eV), while 
for Ei = 9001 eV, we find 1.6 + 0.5% (6 eV). Only the last 
three 4-fold oscillations, albeit small, are statistically signif- 
icant (better than 3cr) and inconsistent with being statistical 
noise. Fourier analysis of the (/f-scan at 3.5 eV (E, = 8997 eV) 
is presented in Fig. [If, where the amplitude is in percent of 
the 0-fold (DC) component. Of the four physical components 
(0-,1-,2-,4-fold), only the 0-, 2-, and 4-fold components are 
above the statistical noise level. The 2-fold component is 
systematically observed in all i/'-scans, but filtered out by the 
symmetrization procedure. 

We can rule out the extrinsic effects of the experimental 
configuration details as the origin of the observed 4-fold os- 
cillations. Extrinsic effects due to sample self-absorption and 
the X-ray beam 'footprint' are expected to be 1- and 2-fold 
respectively. In order to minimize these effects, we aligned 
the scattering surface normal along the i/'-scan rotation axis 
(i.e., (jj -Q in Fig. [T^) and we made sure that the X-ray beam 
spot was contained by the sample surface. After this proce- 
dure, no measurable 1-fold component was observed, but a 
systematic 2-foId component of 2% remained. This 2-foId 
component is observed both within our RIXS signal and with 
a fluorescence monitor placed below the analyzer. Its strong 
dependence on the X-ray beam position within the scattering 
surface confirms that it is a footprint effect. Furthermore, ex- 
trinsic effects should be independent of both the selected in- 
termediate state and the final electronic state that is probed 
and, as a result, be equally prominent for all energy-loss and 
incident-photon energies. This is inconsistent with the vari- 
ability of the 4-fold oscillation amplitude we observe. The 
rotation of the anisotropic momentum-resolution-ellipsoid in 
reciprocal space (see Fig[T}5), might create an artificial 4-fold 
oscillation pattern in the i/r-scans, but a calculation based on 
the ellipsoid's shape and on the anisotropy of the momentum 
dependent inelastic signal gives an upper bound on this effect 
to less than 0.05%, much smaller than the observed oscillation 
amplitudes. 

We can furthermore rule out that the 4-foId oscillations re- 
sult from the resonant nature of RIXS, i.e., from the details 
of the resonantly excited intermediate state. In our experi- 
ment, the intermediate state consists of a li core hole and 
an in-plane Ap photoelectron created by the absorption of the 
incident photon. A local distortion of the lattice could intro- 
duce a 2-fold oscillation by splitting the Ap^ and Apy levels, 
as in the case of the manganites,^^ but the crystal structure 
of Nd2Cu04 is tetragonal and the Ap^- and Apy levels are de- 
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generate. The 4-p states mix non-locally with the Sd^i^yi state, 
which introduces a 4-fold component to the intermediate state. 
On the other hand, two factors make this effect negUgible: 
quadrupole transitions are typically two orders of magnitude 
weaker than dipole transitions and the 3d-4p mixing is weak 
because of small overlap integrals and a large energy sepa- 
ration between the two bands. In addition, an intermediate- 
state effect should be independent of energy-loss value and 
be seen equally in all measured spectra which is inconsistent 
with what we observe. 

We conclude that the 4-fold oscillation is a property of the 
electronic excitations in the final state. That is, as we vary the 
energy-loss value, we probe varying admixtures of final states 
with different symmetries. 

Since the incident-photon polarization is kept within the 
Cu02 planes during the azimuthal rotations, the 4/9-as- 
spectator approximation predicts the RIXS signal to be inde- 
pendent of tff. While this accurately describes approximately 
95% of the inelastic signal, the observed 4-fold oscillations 
with peak-to-peak amplitudes up to 6.4% are evidence of exci- 
tations created by the interaction of the 4p photoelectron with 
the valence system. Because the 4-p photoelectron can trans- 
fer angular momentum, these excitations can be of a differ- 
ent nature than those created by the 1 s core-hole. In indirect 
K-edge RIXS, this result constitutes the first evidence of scat- 
tered photon polarization effects beyond the 4/7-as-spectator 
approximation. 

In Fig. [T^, we extend our analysis beyond the zone cen- 
ter This figure shows the i^-angle dependence of the inelastic 
specti-a at (0 tt) at E, = 8997 eV and 9001 eV. The molecular 
orbital excitation at 5.6 eV measured for E, = 9001 eV has 
a weak iff angle dependence: it is approximately 6% weaker 
for 1^ - 45° than for i/' - 10.7°. There might also be a in- 
dependent feature for E, = 8997 eV at approximately 3.8 eV. 



B. Scattering-geometry dependence 

Scattering-geometry dependence of the RIXS signal can 
have three different origins: sample self-absorption, incident 
and scattered photon polarization effects, and momentum- 
transfer effects. Since all three generally affect the signal si- 
multaneously, it is difficult to separate their contributions. In 
Sec. Ill Al we used a particular scattering geometry that al- 
lows for the polarization degrees of freedom to be varied in- 
dependently of the momentum transfer while minimizing self- 
absorption effects. Alternatively, the scattering geometry de- 
pendence can be investigated by measuring the Brillouin zone 
dependence of the inelastic signal. 

Recently, Kim et al. studied the Brillouin-zone depen- 
dence of charge-transfer excitations at the Cu K-edge of the 
Mott insulator La2Cu04 . Based on remarkable agreement 
between inelastic spectra taken at high symmetry points of 
different Brillouin zones, they concluded that the RIXS sig- 
nal is the same in all Brillouin zones. However, a closer look 
at the experimental spectra shows an interesting energy-loss- 
dependent difference of approximately 10% between spectra 
measured in different Brillouin zones. 
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TABLE I. Normalization factors for six different scattering geome- 
tries. The raw fluorescence monitor (FM) normalization factors are 
compared to different correction methods (see text for details). 



In this Section, we study this subtle effect further by mea- 
suring the scattering-geometry dependence of the related Mott 
insulator Nd2Cu04. Inelastic line scans (energy-gain side sub- 
tracted) taken at the zone center of six different zones are 
shown in Fig. |2^. The proper normalization of each spectrum 
is important as it allows the comparison of relative intensities 
even between widely different scattering geometries. 

In order to separate extrinsic effects from intrinsic features, 
we compare four different normalization techniques. The first 
(FM) consists of using the fluorescence signal as a monitor 
The second (FM-hSA) adds the self-absorption correction de- 
scribed in Appendix |A] The third (SWl) and fourth (SW2) 
both use the fluorescence signal as a monitor and further nor- 
malize the spectra by the integrated inelastic spectral weight 
in a fixed energy-loss range, between 1.4-2.9 eV and 1.9-2.4 
eV respectively. The resulting normalization factors are com- 
piled in Table U where the 'raw' fluorescence monitor factors 
are separated from the different corrections using either self- 
absorption or integrated spectral weight. As expected, (0 L) 
scattering geometries (where w = 0) all have approximately 
the same normalization factor, except for (0 7.1), where the 
X-ray footprint starts to be limited by the scattering surface 
size. The (H L) grazing-incidence normal-emission geome- 
tries (where lj < 0) suffer less self-absorption and have lower 
normalization factors accordingly. Based on the tabulated val- 
ues, the fluorescence monitor provides the largest contribution 
to the normalization 20%), which validates its use as a 
first-order self-absorption correction. However, the additional 
self-absorption correction described in AppendixlAlis not neg- 
ligible 5%) and should be used. 

We analyze the scattering geometry dependence by divid- 
ing each of the six energy spectra into energy-loss bins, each 
bin half the size of the experimental resolution (0.125 eV). 
Within each bin, the scattering-geometry dependence of the 
intensity is fitted to a linear form (y - cri + (T2x), where the 
dependent variable (x) is either L, \Q\^, orcos^(0-w) (see Fig. 
[1^ for the definition of 9 and w). Each fit is performed using 
the four different normalization procedures described above. 

The quality of each linear fit is characterized by a reduced 
chi-square (the total chi-square divided by the number of 
degrees of freedom d). The chi-square values and degrees of 
freedom can be summed across the energy bins to compose a 
collective reduced chi-square - YiiX] 1 2/ (/: bin index) 
which characterizes the fit functions' ability to represent the 
observed Brillouin-zone dependence across the entire data set. 



5 




Energy loss (eV) Energy loss (eV) 



FIG. 2. (a) Two-dimensional zone-center spectra collected for a wide 
range of scattering geometries. Data sets are separated vertically by 
0.3 units. The solid lines through the data are smooth interpolations 
and they are also overlaid at the bottom of the figure to highlight the 
evolution of the scattering geometry data. The data sets are scaled 
using the FM+SA method (see text for details), (b-c) Fits of the 
Brillouin-zone dependence to the linear form y = (T[ +0-2 cos^{9 - ai) 
at 1.65, 3, 4, and 6.5 eV, comparing the FM+SA normalization pro- 
cedure with (b) FM and (c) SW2 normalizations, (d-e) Fitted coeffi- 
cients (Ti and 0-2 as a function of energy-loss, comparing the FM+SA 
normalization procedure with (d) FM and (e) SW2 normalizations. 



The goodness-of-fit indicator (G) is also calculated.^ 

The values of these two indicators, for each combination of 
fit function and normalization procedure, are compiled in Ta- 
ble The cos^(6 - w) linear dependence provides the best 
fit to our data as it robustly yields the lowest and the high- 
est G values, independently of the normalization procedure. 
The L and \Q\^ fits are poorer, no matter what normalization 
procedure is used. 

Examples of the linear cos^(0 - a)) fits are presented in Fig. 
I2J3-C and compare different normalization procedures. The 
energy-loss values chosen for this comparison span the spec- 
tral range of our data and correspond to the thick and gray 
vertical lines in Fig. |2^,d-e. In Fig. |2}i-e, the resulting fit pa- 
rameters (intercept cri and slope 0-2) for each energy bin are 
compared for three of the four different normalization proce- 
dures. 

The choice of normalization procedure affects the extracted 
amount of Brillouin zone dependence. For example, the 
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TABLE II. Reduced chi-square and goodness-of-fit (G) values 
for each combination of fit function and normalization procedure. 



FM+SA coiTection in Fig. adequately includes the ef- 
fects of scattered-photon polarization: the self-absorption is 
larger for normal emission (in-plane scattered photon polar- 
ization, cos^(6 -(L)) ^ Q) than for grazing emission (mix of in- 
plane and out-of-plane photon polarization, cos^(6 - w) 1). 
This difference is apparent form the polarization dependent 
X-ray absorption curves in Fig. H):. 

The SW2 normalization yields the lowest ™d largest 
G which suggests that the 1.9-2 A eV energy-loss region is 
Brillouin-zone independent. On the other hand, this nor- 
malization produces unphysically high values of 0-2 (Ref.— ) 
which suggests that the 1-3 eV energy-loss region in 0-2 
is artificially reduced by a source of error beyond the self- 
absorption correction. 

While a priori not unphysical,^' the negative offset in 1x2 
below 3 eV (using FM+SA normalization) is probably an ar- 
tifact the elastic-line subtraction. In our subtraction of the 
energy-gain side, we assume that the elastic line is symmet- 
ric. While this approximation is in principle valid, weak 
anisotropic elastic signal could 'leak' into our RJXS spectra 
(due to the non-zero energy resolution) and introduce an ar- 
tificial Brillouin-zone dependence of the signal at low energy 
loss. Errors in the fitted elastic line position can also introduce 
a weak Brillouin zone dependence. While this anisotropy is 
limited to low energy loss, an eiTor in fitted position would 
be proportional to the RJXS spectrum's energy-loss slope and 
create artifacts at both high and low energy-loss. We em- 
phasize that whereas the broad negative offset observed in 
(72 below 3 eV (within FM+SA) could be explained by a 
Brillouin-zone-varying asymmetry in the elastic line, neither 
of the above effects can create sharp features like the one ob- 
served at 1.65 eV. 

Finally, a slight crystal misalignment could account for the 
overall larger intensity measured between 3-5 e V at (0 11 . 1 ), 
since the zone center spectrum at 8997 eV is a local minimum 
of inelastic intensity. On the other hand, the different spectral 
shape observed for (2 11) and (3 12) around 4 eV is not 
explicable by a crystal misalignment and could be evidence 
of a Brillouin-zone dependence not captured by the analysis 
presented in Fig. |2] 

While the many sources of eiTor make the quantitative com- 
parison of the spectra difficult, the quantities cri and 0-2 ex- 
hibit robust features. The Brillouin-zone independent part of 
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the RIXS spectrum (tri) consists of a broad feature centered at 
4.5 e V, a shoulder at 2 eV (and possibly another at 1 .4 eV), and 
a (linearly-extrapolated) onset of 0.8 eV. On the other hand, 
the Brillouin zone dependent part (0-2) consists of a broad 
feature, centered around 5 eV with an onset around 2.5 eV, 
and its most interesting feature is a resolution-limited peak at 
1.65 eV. 

Because the spectra cannot be scaled to collapse onto one 
common curve, no matter what normalization is used, the 
Brillouin-zone dependent part (1x2) cannot be spurious. The 
energy-loss dependence also rules out the resonant cross- 
section as the origin of the Brillouin-zone dependence, leav- 
ing only the properties of the measured electronic excitations 
to explain the effect. 

For the employed scattering geometry, the 4/9-as- 
spectator approximation predicts no photon-polarization- 
based Brillouin-zone dependence of the inelastic signal. In 
contrast, the observed Brillouin-zone dependence is best fit by 
the cos^(6 - a>) form which implies that the effect is photon- 
polarization-based and not momentum-based. As such, this 
observation implies that the 4p photoelectron interacts with 
the valence system during the RIXS process. 

We note that the integration of the FM-hSA normalized 
0-2 spectral weight above 2.5 eV-'^ sums to 15 + 2% of the 
integrated cri spectral weight. This departure from the 4p- 
as-spectator approximation is larger than the i/r-scan peak-to- 
peak amplitude variations (~ 6%), but it is of the same order 
of magnitude, which suggests a similar mechanism for both 
effects. 



III. DISCUSSION 

The analysis of the incident- and scattered-photon polar- 
ization dependence of the cross section is key to the study of 
the symmetry of electronic excitations in Raman scattering. 
Within (direct) soft RIXS, this has been used to distinguish 
the Zhang-Rice singlet (ZRS) excitation from local d d 
excitations.^'^'* On the other hand, for (indirect) hard RIXS, 
it is unknown what excitation symmetries are measured and 
what their relative strengths are. It is furthermore unknown if 
RIXS obeys selection rules linking the incident- and scattered- 
photon polarizations and the underlying excitation symme- 
tries. 

While theoretical treatments have made assumptions about 
what types of excitations are measurable, they have not dis- 
cussed selection rules explicitly. Treatments using joint- 
density-of-states-type cross sections limit the scattering from 
the core-hole to interband transitions between bands of the 
same point-group symmetry at g = 0.^^'^^ Furthermore, cal- 
culated RIXS spectra using one-band Hubbard models are au- 
tomatically limited to states of - local symmetry (the 
ZRS combination of oxygen orbitals is an - combina- 
tion of the O Pa- orbitals), so that the symmetry of the charge- 
transfer excitations is limited to Ai^ at g - 0?^ On the other 
hand, non-Ai^ transitions have been suggested to explain new 
features in measured RIXS spectra: charge-transfers to non- 
bonding bands"* and local d ^ d excitations^^ are both exam- 



ples of such non-A]^ transitions. 

While it remains unclear if the polarization-based Raman 
selection rules can be applied to interpret indirect RIXS spec- 
tra at the Brillouin zone center (zero reduced q), we test their 
applicability by comparing their predictions with our data. 
From the definition of the symmetry channels allowed by the 
tetragonal (D4/,) crystal structure of Nd2Cu04 (Ai^, A2g, Big, 
B2g, and Eg) and the evolution of the incident- and scattered- 
photon polarizations with the azimuthal angle (tfr) and the an- 
gular difference (6-(d), we can write the photon-polarization- 
based Raman selection rules, at Q = 0, as a function of the 
scattering power within each allowed symmetry channel: 

/inel. ^ [o-(Aig) + 0-(A2g) + O-(Big) + 0-(B2g)] 

+ [cr(Eg) - cr(Big) - cr(A2g)] cos^(0 - u) 

+ [o-{Big) - o-(B2g)] cos2(2iA) cos\d - cj) (2) 

Since the scattered-photon polarization is not analyzed, we 
must include both the cr and n channels (polarization perpen- 
dicular and parallel to the scattering plane, respectively), as 
shown in Fig. [T^,b. Note that these rules inform us about the 
polarization dependence of these symmetry channels, but not 
about their relative size. 

At fixed 6 and for lj - Q, Eq. |2] can be rewritten as 
(To -I- 0-4 cos^(2(/') which is precisely the functional form of the 
4-fold oscillations presented in Sec. Ill Al With a peak-to-peak 
4-fold amplitude equal to |^cr(Bi^) - cr(B29)] cos^(9-a)), these 
azimuthal oscillations can be interpreted in terms of the Big 
and B2g symmetry channels. While we cannot determine the 
Big and B2g amplitudes independently, our data suggest the 
presence of Bi,j-type electronic excitations at 3.5 and 5 eV en- 
ergy loss for Ei = 8997 eV and at 6 eV for = 9001 eV. Cor- 
recting the 4-fold peak-to-peak amplitude to account for the 
cos^(6' - ai) dependence, we obtain an adjusted Big - B2g am- 
plitude of +8.7% (0-4/0-Q) between 3.5 - 5eV for E,=8997 eV. 
These adjusted data are shown in Fig. [3] 

For i// - 0, Eq. |2] reduces to a cross section of the form 
(Ti + (72 cos^{9 - cj), where cri - Aig + A2g + Big + B2g and 
0-2 - Eg - A2g - B2g, which is precisely the function that best 
fits the Brillouin-zone dependence presented in Sec. IIIBI This 
suggests that the decrease in spectral weight toward backscat- 
tering observed here for Nd2Cu04 and by Kim et al. for 
La2Cu04 can be interpreted in terms of Eg excitations at high 
energy-loss. These symmetry assignments for cri and 0-2 are 
reproduced in Fig. [3] 

In tetragonal symmetry, the incident- and scattered-photon 
polarizations of an Eg excitation correspond to 4p crystal-field 
eigenstates with different energies (e.g., x— >z or z^x), unlike 
the Aig, A2g, Big, B2g symmetry channels where the incident 
and outgoing 4p crystal-field eigenstates are degenerate (and 
planar). This creates a unique resonance profile '^i'^'^^i^^ for 
Eg excitations which an increase or decrease of intensity com- 
pared to the other symmetry channels. 

For example, in the scattering geometry used here, at the 
incident energy E, = 8997 eV, Eg and Big excitations have 
the same incident-photon polarization resonance but their 
scattered-photon polarization resonances differ. For Big ex- 
citations, the scattered-photon polarization is e/ = x and the 
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FIG. 3. Symmetry interpretation of the photon-polarization depen- 
dent and independent RIXS spectral weight for tetragonal Nd2Cu04. 
While composed of many different symmetry channels, the blue 
spectrum is most likely dominated by the Aig symmetry channel. See 
text for details. 



final energy resonance is at - 8997 eV whereas for Eg 
excitations, the scattered-photon polarization is e/ = z and 
there are scattered-photon resonances at Ef - 8985, 8993, 
and 8997 eV according to the XAS data in Fig. |4] Around 
4 eV and 12 eV energy-loss, the scattered photon is tuned to 
an out-of -plane resonance and Eg excitations will be enhanced 
compared to Big excitations. In this case, assuming that both 
symmetry channels have a comparable density of states, we 
estimate that Eg excitations are comparatively enhanced by 
a factor of 2-5. If we correct for this enhancement, the ob- 
served 15% Brillouin-zone dependence becomes a 3-8% ef- 
fect, more comparable to the observed 8.7% azimuthal-angle 
dependence. 

For the scattering geometry used here, the core-hole is ex- 
pected to create the majority of the inelastic signal because 
it has the strongest effect on the valence system. We suggest 
that the Brillouin-zone-independent and i/'-angle independent 
spectral weight is created by core-hole scattering and consists 
ofAig symmetry. We note, however, that this contribution can 
be reduced to zero in certain scattering geometries because its 
cross-section follows the 4/5-as-spectator approximation dis- 
cussed in Sec. I] 

While we have shown that the Raman selection rules can 
accurately model the observed photon-polarization effects and 
suggest symmetry assignments for different inelastic features, 
these assignments should be supported by a theoretical un- 
derstanding of the excitations. The molecular orbital (MO) 
excitation at the poorly-screened state is understood to be of 
A[g symmetry, which is consistent with its weak i^-angle de- 
pendence at (0 n) (see Fig. [1^). While the MO excitation does 
show a weak 4-fold oscillation at the zone center, a compari- 
son of the (0 n) and (0 0) spectra suggests that the 4-fold oscil- 
lation might instead be a property of a momentum-dependent 
shoulder that disappears away from (0 0). 

In principle, d ^ d excitations provide a testing ground 
for the validity of conventional Raman selection rules within 
RIXS, and such excitations have been well studied with soft 



X-ray RIXS at the Cu L-l-iOiill and M-edges," with opti- 
cal absorption,"*^ with large-shift Raman scattering,'*^ and with 
different theoretical methods.^^^^ 

Although for Nd2Cu04 only the A2g (d^2_y2 -> dxy) Cu 
crystal-field excitation has been observed around 1 .4 eV,"*'' the 
crystal-field excitations of tetragonal Sr2Cu02Cl2 have been 
extensively studied and their energies should be sim- 

ilar to those of Nd2Cu04. For the latter, Raman selection 
rules suggest that the resolution-limited feature we observe 
at 1 .65 eV (Fig. ^ is of Eg symmetry. The Eg crystal-field ex- 
citation (d^i^yi djc,) in Sr2Cu02Cl2 has an energy of 1 .7 e V 
which, supports this symmetry assignment and the applicabil- 
ity of the Raman selection rules. 

A complete determination of the d d excitations in 
tetragonal Nd2Cu04 should be possible with the experimen- 
tal methods described in this paper, i.e., without analyzing the 
scattered-photon polarization, since only three types of d d 
excitations are possible (the B2g symmetry change does not 
exist within the Cu 3d orbitals). For example, with better 
energy resolution, the A2g excitation should be observable at 

1 .4 eV with a Brillouin-zone dependence study, whereas the 
Big {dxi_yi — > d^,2_i.i) excitation should be observable with a 
i/'-angle dependence study. 

Charge-transfers from the non-bonding oxygen bands to 
the upper Hubbard band, which have been used to explain 
the multiplet of inelastic features between 2 and 6 eV^ 
are expected to be of non-Ai^ symmetry. As studied with 
ARPESr^^ O 2p non-bonding bands have approximately 

1 .5 eV more binding energy than the ZRS band. Accordingly, 
the RIXS charge transfer to such bands should start at an en- 
ergy 1.5 eV above the overall onset of RIXS excitations. In 
Nd2Cu04, the onset of excitations is approximately 0.8 eV. 
On the other hand, the 4-fold azimuthal-scan oscillations start 
between 2 and 3.5 eV, and the broad feature in the Brillouin- 
zone dependent part starts at around 2.5 eV. Both have onsets 
approximately 1.5 eV above 0.8 eV and could be interpreted 
as Big and Eg charge-transfers to oxygen non-bonding bands, 
respectively. 

In their EELS study of Sr2Cu02Cl2, Moskvin et al.^ iden- 
tify many charge-transfer excitations, three of which have Ai^ 
symmetry and should be RIXS active. The two strongest exci- 
tations are around 8 eV, but do not have obvious RIXS equiv- 
alents. The third Aig excitation at 2 eV, just above the optical 
gap, is identified as the Zhang-Rice singlet. This may cor- 
respond to the 2 eV feature in RIXS, which is apparent as 
a weak shoulder in Fig. |2]l-e and is seen as a clear peak in 
La2Cu04 and Sr2Cu02Cl2. This suggests that the 2 eV fea- 
ture in RIXS is of Aig symmetry. An exact diagonalization 
calculation Tohyama of charge-transfer excitations within a 
one-band Hubbard model of the Cu02 plane agree quite well 
with our data and the Aig symmetry assignment of the 2 eV 
feature. The calculation shows a peak at 2 eV that is exclu- 
sively of Aig symmetry, and a continuum of excitations be- 
tween 2 - 3.5 eV that is of predominantly Bi^ character, in 
good agreement with the symmetry interpretation of our data. 
The relative intensity of these two features disagrees with our 
observations, but quantitative agreement is not expected, since 
the resonant and non-resonant cross-sections are very differ- 
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ent. 

In order to better understand the symmetry selectivity of 
the RIXS cross section, it is worthwhile to compare it to well- 
established probes. On the other hand, symmetry-selectivity 
differs from probe to probe, which renders direct compar- 
isons hazardous. For example, in crystals with an inversion 
center, optical conductivity measures a current-current corre- 
lation function only sensitive to ungerade excitations, while 
electron-energy-loss-spectroscopy (EELS) and IXS measure 
a density-density correlation function that is sensitive to both 
ungerade and gemde excitations. Furthermore, the symmetry- 
sensitivity of these probes has a strong and non-periodic Q- 
dependence^''^^ in contrast to the nearly Brillouin-zone inde- 
pendent RIXS spectra. 



IV. CONCLUSION 

By studying the photon-polarization dependence of indi- 
rect RIXS at the Cu A'-edge of the tetragonal Mott insulator 
Nd2Cu04, we uncover anomalous excitations of two differ- 
ent types, which is evidence of scattered-photon polarization 
effects. While the majority (80%) of the inelastic signal is 
describable by the 4/?-as-spectator approximation, the anoma- 
lous remainder constitutes the first evidence of RIXS excita- 
tions created by the interaction of the Ap photoelectron with 
the valence system. 

The successful modeling of the observed azimuthal-scan 4- 
fold patterns and of the Brillouin-zone dependence by photon- 
polarization-based Raman selection rules suggests that these 
rules can be used to interpret zone-center RIXS spectra. Using 
these rules, we tentatively assign the sharp peak at 1 .65 e V in 
the Brillouin-zone-dependent spectral weight to an Eg d 
d excitation and the broad features above 2.5 eV in both 
the Brillouin-zone-dependent and azimuthal-angle-dependent 
spectral weight to Eg and Big charge-transfers to non-bonding 
oxygen bands. 

Establishing photon-polarization-based methods to charac- 
terize the electronic excitations' symmetry is a pivotal chal- 
lenge for RIXS. Such methods should facilitate the interpre- 
tation of experimental spectra and help provide a better un- 
derstanding of the underlying physics of the cuprates and 
other transition metal oxides. While a complete scattered- 
photon polarization analysis is currently prohibited by low 
count rates, ongoing instrumentation development and future 
increases in photon flux should soon make the full polariza- 
tion analysis of the inelastic signal possible. 

We would like to acknowledge valuable conversations with 
J. van den Brink, T.R Devereaux, and K. Ishii. This work 
was supported by the DOE under Contract No. DE-AC02- 
76SF00515 and by the NSF under Grant No. DMR-0705086. 



Appendix A: RIXS Normalization with a fluorescence monitor 

Aside from effects intrinsic to the cross section, the RIXS 
scattering intensity is also modulated by the sample self- 
absorption, an extrinsic effect that depends on the scattering 



geometry. 

In X-ray scattering, the absorption processes which deter- 
mine the X-ray absorption length of a crystal are dominated 
by Auger emission of electrons and core-level fluorescence 
lines. Together, these processes' cross-sections dwarf the 
RIXS cross-section. As a result, when the X-ray beam travels 
through the sample, before and after the RIXS event, its inten- 
sity is attenuated. This attenuation has a strong dependence on 
the scattering geometry and determines the illuminated sam- 
ple volume. The RIXS signal is proportional to the number of 
Cu atoms resonantly excited in this volume. This attenuation 
is referred to as self-absorption, and correcting for it is com- 
mon practice for other spectroscopies scattering probes (for 
example, in EELSjS in EXAFS,21 and in direct RIXS^^ for 
example) but has not yet become standard for indirect RIXS. 

To calculate this effect, we must know the X-ray absorption 
length and the scattering geometry. In reflection geometry, 
the RIXS intensity given by Eq. lAll It includes the incident 
beam intensity /o, the intrinsic RIXS scattering amplitude per 
Cu atom F{Ei, AE), the incident beam cross-sectional area Ba, 
the density of Cu atoms pcu, the absorption coefficient for the 
incident (scattered) photons yu,- (jjf) with polarization e, (Cf) 
and energy £, (Ef), as well as the scattering angles 0, and 
Of given relative to the sample surface, as shown in Fig. I^^. 
The maximum scattering intensity is obtained for a photon 
with grazing incidence which is emitted perpendicular to the 
sample surface (an example is shown in Fig. In turn, 

the minimum intensity is observed for normal incidence and 
grazing emission angle. The reduction in RIXS signal from 
having a footprint (the beam spot on the sample) larger than 
the scattering surface is not included but is an important effect 
at grazing incidence. 

/ = (kpcuBa) S [Ei, ei,Ef, ef) F (£,■, - Ef) 

Me„£,)+Mf/,£/)i^ ^ ^ 

In an ideal experiment, this formula could be used directly 
to normalize the RIXS spectra taken in different scattering ge- 
ometries. In practice though, the footprint is highly depen- 
dent on the scattering surface and can affect the RIXS signal 
in non-trivial ways. To counteract this difficulty, we use a 
fluorescence detector put at a known position close to the an- 
alyzer crystal, as presented in Fig. H^. After measuring the 
position of the sample surface relative to the crystal axes and 
tuning the fluorescence detector to an emission line (for ex- 
ample Cu Kai), we can calculate the scattering geometry de- 
pendence of both the fluorescence monitor (FM) signal (S fm) 
and of the RIXS signal (S rjxs )■ In order to correct for the self- 
absorption effect, which changes the intensity of the RIXS 
signal based on scattering geometry, we normalize with the 
following ratio: 



I RIXS S FM 

X 



' RIXS 



(A2) 



The measured fluorescence and RIXS intensities are Ifm and 
Irixs ■ Simply dividing the RIXS signal by the FM signal is 
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FIG. 4. (a) Position of the fluorescence monitor relative to tiie an- 
alyzer crystal (b) Self-absorption parameters in reflection geometry, 
(c) X-ray absorption spectroscopy (XAS) data (by partial fluores- 
cence yield of the Cu Ka\ emission line) in the Cu A'-edge region. 
The black lines correspond to the calculated X-ray absorption coeffi- 
cients above and below the Cu A'-edge for Nd2Cu04. The measured 
XAS curves are scaled to follow the published curves by Tranquada 
et al.S' . (d) Nd2Cu04 absorption curve at the Cu A'-edge measured 
with the fluorescence monitor as described in the text. The black 
arrows indicate incident energies of 8997 eV and 9001 eV. 



a first-order correction that can partially account for changes 
in the X-ray footprint on the sample. However, the exact cor- 
rection must include variations in X-ray self-absorption based 
on the scattering geometry, the photon's polarization and en- 
ergy all of which require calculating the self-absorption fac- 
tors Stm and Srixs- Note that the S fm calculation includes 
the difference in location between the spectrometer's analyzer 
crystal and the fluorescence monitor, as shown in Fig. |4^. In 
practice, the further away the FM is placed from the analyzer 
crystal, the more difficult it will be to accurately account for 
positional differences in the calculation of the self-absorption 
correction. This can become a large source of systematic error 
in the normalization procedure, although this is not the case in 
the present experiment. 

The fluorescence monitor signal described above is col- 
lected with a solid-state detector (Amptek). The signal is 
fed into a scaler and only photons within a fixed energy 
range are counted, in our case, in a 1 keV range around the 
strong Cu Kai emission line at 8.1 keV. An absorption curve 
for Nd2Cu04 with in-plane incident-photon polarization col- 
lected with this fluorescence monitor is presented in Fig. |4}i. 
The non-zero signal below the absorption edge is parasitic but 
is not the result of dark current within the detector Instead, 
it is composed of two signals that leak into the energy inte- 
gration window: elastic scattering at 9 keV and neodymium 



L-edge fluorescence at 7 keV. Within our theoretical calcula- 
tion, this parasitic component is tentatively compensated for 
by modeling it as 70% elastic signal and 30% neodymium 
fluorescence. Since the parasitic component does not vary 
congruently with the measured emission line, it constitutes a 
source of error and should in general be minimized. 

In order to calculate the polarization dependent X-ray ab- 
sorption coefficients (fi) in the Cu K-edge region, we collected 
X-ray absorption data by partial fluorescence yield at the Cu 
Ka\ emission line with an X-ray spectrometer^^ (1.2 eV reso- 
lution) The effect of self-absorption, typically important close 
to absorption edges, was calculated as described by Carboni 
et al. , but was not significant here. In order to complete the 
construction of the absorption coefficients, the curves were 
then scaled to match the tabulated values for X-ray absorption 
in Nd2Cu04 above and below the Cu /T-edge and following 
Ref.^^. The results are shown in Fig. |4j;. While both the nor- 
malization procedure above the edge and the self-absorption 
correction of the fluorescence signal are sources of error, a 
conservative estimation of their effects is included in the nor- 
malization factor's error but, in our case, they are negligible. 

This type of correction includes absorption effects due to 
changes of the scattered-photon energy that simply dividing 
by the fluorescence signal cannot account for For example, 
based on Fig. H};, an incident in-plane photon with an energy 
of 8997 eV would see more absorption on its way into the 
sample than the scattered photon on its way out of the sample 
with 10 eV energy-loss (Ej - 8987 eV). Within an inelastic 
spectrum taken in these conditions and corrected by the fluo- 
rescence signal alone, the high-energy-loss response would be 
artificially larger than that at low energy, based solely on this 
change in the XAS coefficient. 

Finally, the scattered photon can in principle have cr or n 
polarization which affects the calculated self-absorption cor- 
rection. In vertical scattering geometry, at the in-plane Ap 
resonance, and within the bounds of the 4/?-as-spectator ap- 
proximation, the scattered photon must be cr polarized and we 
calculate the self-absorption correction based on this prescrip- 
tion. As seen in Sec. IIIBI there is evidence that this approxi- 
mation only describes part of the RIXS signal, with the com- 
plementary contribution in the n polarization channel. This 
latter polarization channel, when isolated, can be normalized 
by modifying the calculated scattered photon X-ray absorp- 
tion coefficient (jij) accordingly and multiplying the cr — > tt 
scattering intensity by the ratio of the cr ^ n and cr — > cr 
calculated self-absorption corrections: 



■cr— >;r 
RIXS 



(A3) 



On the other hand, within the measured spectral range, the 
maximum correction is only 4% of the Brillouin-zone depen- 
dent signal and does not noticeably affect the normalization 
outlined above. 
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